One-loop renormalization group study of boson-fermion mixtures 
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Abstract 



A weakly interacting boson-fermion mixture model was investigated using Wisonian renormalization group analysis. 
This model includes one boson-boson interaction term and one boson-fermion interaction term. The scaling dimensions 
of the two interaction coupling constants were calculated as 2 - Z) at tree level and the Gell-Mann-Low equations were 
derived at one-loop level. We find that in the Gell-Mann-Low equations the contributions from the fermion loops go to 
zero as the length scale approaches infinity. After ignoring the fermion loop contributions two fixed points were found 
in 3 dimensional case. One is the Gaussian fixed point and the other one is Wilson-Fisher fixed point. We find that the 
boson-fermion interaction decouples at the Wilson-Fisher fixed point. We also observe that under RG transformation the 
boson-fermion interaction coupling constant runs to negative infinity with a small negative initial value, which indicates a 
boson-fermion pairing instability. Furthermore, the possibility of emergent supersymmetry in this model was discussed. 

Introduction 



the interspecies attraction in boson-fermion mixtures was 
observed experimentally |9|. Also, the formation of het- 

eronuclear Feshbach molecules has been observed in a boson- 

in "^He in 1995111, the field of degenerate quantum gases ^ • • . csini. j 4()t/^ . • or-. 

'^ ^ '^ rermion mixture or Kb and K atomic vapors in a 3D op- 
has become one of the most active areas of physics. Of ,. , , ,.,. JTrn j • 4.- i j- i 4_ irm 

^ •' tical lattice BlOl and in an optical dipole trap lllll . 



Since the first observation of Bose-Einstein condensation 

4 



particular interest is the realization of boson-fermion mix- 
tures of atom gases. They may show very different behavior 
from pure fermion or pure boson gases. Various theoret- 
ical researches have been proposed. For instance, forma- 
tion of stable strongly correlated boson-fermion pairs(2], 
instability of the mixture when there is an attraction be- 
tween bosons and fermions|l3]|4l, interspecies interactions 
induced attraction among bosons|l5]|6l and emergent super- 
symmetry (SUSY) from mixtures of cold Bose and Fermi 
atoms||2l[8l. Recent developments in atomic experiments 
have made it possible to realize boson-fermion mixed gases 



In the present work we give a renormalization group 
analysis on a boson-fermion mixture model at finite tem- 
perature. Wilsonian renormalization group approach [ll4] 
[TSl is a popular method to study various condensed mat- 
ter problems. This technique has been applied to a homo- 
geneous Bose gas by several authors llT6l [TT] [TSl . How- 
ever, it was recognized in 1990s that the standard Wilson's 
momentum-shell approach must be modified for systems in- 
volving Fermi surface lfT9l l20l UA] since in such a system we 
renormalize not towards a single point, the origin, but to- 
wards the Fermi surface. Renormalization only reduces the 



in the laboratory. Collapse of the atomic cloud induced by ,. . i ^ .u t- • t u-i .u . .■ i 

■> t- J dimension normal to the Fermi surface while the tangential 



part survives ll22l . Besides the applications of renormaliza- 
tion group in pure-boson and pure-fermion systems a RG 
formalism for mixed boson-fermion systems were also dis- 
cussed by several authors ||23] l24l l25l l26l l27l . 

In this context one is dealing with dilute, weakly in- 
teracting systems. This allows to effectively express the 
quantities of interest in terms of a single parameter char- 
acterizing the particle interaction. Our boson-fermion mix- 
ture model includes two important interaction parameters 
gi and g2 which denote short-range boson-boson interaction 
and boson-fermion interaction respectively. The renormal- 
ization group analysis shows that the scaUng dimensions of 
gi and g2 are both 2 - D at tree level, where D is the di- 
mension of the system. Hence, gi and g2 are both marginal 
when D - 2 and irrelevant when D > 3. At one-loop level 
we derived the Gell-Mann-Low equations and found that in 
these equations the contributions from the fermion loops go 
exponentially to zero as /' — > oo compared with the contri- 
butions of the boson loops. After we ignore the contribu- 
tions of fermion loops, two fixed points are found in 3 di- 
mensional case. One is the trivial Gaussian fixed point, the 
other one is the Wilson-Fisher fixed point. At the Wilson- 
Fisher fixed point the parameter g2 goes to zero. This im- 
plies that at one-loop level the boson-fermion interaction of 
this model decouples at the critical temperature. We also 
find that the the boson-fermion interaction coupling con- 
stant with a small negative initial value runs to negative in- 
finity under the renormalization transformation. This could 
indicate a boson-fermion pairing instability. 

In the low-energy limit of a nonsupersymmetric con- 
densed matter system supersymmetry(SUSY) can dynami- 
cally emerge at a critical point lfT2l . For our model if the 
chemical potentials of boson and fermion are equal and the 
two coupling constants are identical the Hamiltonian is in- 
variant under supergroup t/( 1 1 1 ) lfT3ll . We use RG to explore 
if there is such a SUSY fixed point. It turns out in the weak 
interaction limit this model doesn't exhibit a SUSY fixed 
point. 



2 The Model 

The model we concerned with includes one boson field (p 
and one spinless fermion field \j/. The grand partition func- 
tion can be expressed as a functional integral. 
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where 
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We work in D-dimensional space, where the fields depend 
on spatial coordinates x - (xi,X2, ...xd) and the imaginary 
time T. In this paper we consider the cases of D > 2. The 
coupling constants for the short-range boson-boson interac- 
tion and boson-fermion interaction are denoted by g\ and 

gi- 

In order to discuss the scaling of the momentum we ex- 
pand the fields in Fourier modes though 



^«-) = ^Z/^/^^'-'> 
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j\ui-i-uylT) 
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where w* = 2m and w/, - i£21i£L ai-e the Matsubara fre- 
quencies for boson and fermion respectively and/? - llksT. 
kg denotes Boltzmann' constant. Then we can rewrite the 
action in momentum space, 

S[b\bJ,f] 

— -— /(^, u)i){-iJ„ +eK- lif)fii, U) 
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we call [A] the scaling dimension of A. In this manner the 
scaling dimension of momentum q is 



[q\ = 1. 



(8) 



+gi- 
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y, r d'^K4 d'^K2 d'^q-i d'^qi 



f(K4, aji)f(K2, <)fe*(?'3, ^^3)^(^1, w^,) 



■6 (K4 +43-K2-qi)- 5„/^+^<;^,„/^+„j_ 



(5) 



In above equation eg = q^/2mh and €k = K^/2mf are ki- 
netic energies for boson and fermion respectively. 

3 Renormalization Group Analysis 
3.1 Tree Level Scaling 

We follow the Wilson's momentum-shell approach. The 
renormalization group transformation involves three steps: 
(i) integrating out all momenta between A/s and A, for 
tree level analysis just discarding the part of the action in 
this momentum-shell; (ii) rescaling frequencies and the mo- 
menta as (co, k) — > (s^"^a), sk) so that the cutoff in k is once 
again at ±A; and finally (iii) rescaling fields (p — > s^'^^ip to 
keep the free-field action S invariant. 

First Let's think about the quadratic term of the boson 
field. After we integrate out a thin momentum shell of high 
energy mode the limit of q (which is the radial coordinate 
of the momentum space) changes from [0, A] to [0,A/s], 
where s ^ 1 . In order to compare the action with the origi- 
nal one we need to rescale the radial coordinate as 



q = sq. 



(6) 



Hence, the cutoff in q is back again at A. Here we give a 
definition to the scaling dimension. If a quantity scales as 



A' = s'^U, 



(7) 



Then the scaling dimensions of the boson field, the energy 
and the chemical potential can easily be derived from the 
quadratic part of the boson action. Following the first two 
steps of the Wilson's renormalization group transformation, 
the quadratic term of the boson action becomes 

(9) 
To make it invariant under the scaling transformation we 
define the scaling dimension of the boson energy as 

[w*] = 2 (10) 

and the scaling dimension of the boson field as 

D + 2 
[b]^ — . (II) 



Now we turn to the fermion case. The quadratic part of 
the fermion action is given by 



^^Z 



/ ^^^^^' co{,)(-icai+eK-Hf)f(K, co{). (12) 



In contrast to the boson case, low-energy modes of fermions 
live near the Fermi surface. In order to preserve the Fermi 
surface under scaling we can't simply scale the momen- 
tum as we did in the bosonic case. We renormalize not 
towards a single point, the orgin, but towards a surface. 
To make progress we define a lower-case momentum k = 
\K\ - Kf, which corresponds to the low energy mode of 
fermions. Then it is the momentum k but not momentum 
1^1 that scales under the renormalization group transforma- 
tion. Since 



£K-f^f^ 



2m 



-Sfif - Vf(\K\-Kf)-diJ.f = Vfk-d/if, 

(13) 



the quadratic part of the action can be approximated as 



/ 
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jjf{K, u>i){-iu>i + £k- f^f)f(K, w{) 



^n^'K^-' f dkf(k,aji)(-iaji + VFk-6i^f)f(k,U), 

J-A 



(14) 



where v^ is the Fermi velocity and (5/// = f^f(T) - /i/(0) can 
be considered as the chemical potential of the low-energy 
modes of fermions. Following the first two steps of the 
renormalization group transformation this part becomes 

/-A _ , 

n'^K^-^ s-^''^dkf(k',co'/)(-is-^""^co'/ + VFS-^''^k' 

J-A 

-s-^'^f^6i,'f)f(k',co'/). (15) 

In order to analyze fermions and bosons in one model it is 
reasonable to scale the energies of fermion and boson the 
same way, that is 



term. After we throw away the high energy momentum 
shell, the interaction becomes 

■0{Kis - \4A\)e{Ns - m)e{Ns - \q2\)e{Kis - m), 

(20) 



where we implement function to generate constraints on 
the momentum space instead of cutoffs in the limits of inte- 
gration, which gives a more explicit description in the scal- 
ing analysis. We eliminate one momentum variable qn using 
the delta function 6'^(q4+q3-q2-qi)- The above interaction 
term can be written as 



;i 



(2nr 



[w{] = [c^n] = 2. 






(16) 
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According to Eq.(15) the scaling dimension of the low en- 
ergy fermion momentum k is the same as the fermion en- 
ergy, 



n^-' VT'^' "'«2-' 

1+^2 -^3 1) 

■e{kis - m)e{Kis - \42\)0{Ns - m)- 



■'(J'+(J'M"+(J, 



p e(A/s - \4 



(21) 



[k] = [w{] = 2 



(17) 



To take the Eq.(15) back to the original form Eq.(14) we 
have to rescale the fermion field as 



/' - s-^'^f. 



(18) 



Then the scaling dimensions of the fermionic fields is 



[/] = -[k] 



(19) 



So far we have gained the scaling dimensions of mo- 
menta, energies and fields of both boson and fermion. Now 
we are ready to calculate the scaling dimensions of the in- 
teraction coupling constant ^i and g2- The renormaliza- 
tion group transformation of the two-body interaction terms 
shows more subtleties, especially for the boson-fermion in- 
teraction term. First we study the pure boson interaction 



When the momentum qi are scaled as ^J = sqi, the 9 func- 
tions transform as 

6{Kis - m = 0(A - sm 

= 6i(A-|^|), for i=l,2,3, (22) 

and 

e{Kls - \qx +q2- qi\) = 9{K - s\qi + q2 - q^l) 

= 0(A-|^;+^2-?3l)-(23) 

All the 6 functions transform back to the original forms. 
Then we can scale the pure boson interaction term as 



2-Di 



8^ (2nf y r(flf^^ 

2 /?' ^ J m (27Tr 

b'*(4'^ +4',- ^3, <)fe'*(^3, <)fe'(^2, <)fc'(^i, <) 

■e(A - \4',\)0(A - |^2l)0(A - \qi 1). (24) 



Notice that fi scales as the inverse of energy, therefore its but also of Kf. It's easy to check that 0(h - \k4\) doesn't go 
scaling dimension is back to the original one after the RG transformation. 



m 



-2. 



(25) 



0(A - |^4l) 



In order to transform the Eq.(24) back to its original form 
Eq.(20) we define 

(26) 



^^-^^1, 



e(K - (\K2 + 4i- u - Kf)^ 

^ e^KIs - (1^2 + qi- U - ^f)) 
■ e(K-s{\K2 + 4x-qi\-KF)^ 



[K-{\sK2 + q\ -q'^\ - sKp 



)) (31) 



then the scaUng dimension of gi is 
[^i]=2-D. 



As discussed by R. Shankar 112111221 the renormalization 
group transformation of a system involving fermions must 
be treated carelly. Much of the new physics stems from 
measure for quartic interactions involving fermions. The 
boson-fermion interaction term in our model is 

(27r)^ Y r d"K4 d'^K2 d'^qi d'^qi 
^^' B ^ J (2^(2^(2^(2^ 

■d°{K4 +4i-K2- ^l)^„,/^+^i;^,„/^+„i;^ 

■0(A - M)6(K - |^3|)0(A - \k2\)6(K - 1^1 1). (28) 



How can we say what the new coupling constant is if the 
(27) integration measure doesn't go back to its old form? To 
solve this problem we approximate 1^4! as 



1^:4! = \K2 + qi- qi\ - Kf 

= \K2 + 4\-Kf 
^Kf(\K2+4/Kf\-1) 
= Kf(\A\ - 1). 



(32) 



where q - q\ - qj,, |A| = \K2 + q/KF\ and K2 is the unit 
vector of K2. With this approximation the transformation of 
function is written as 

6»(A - \k4\) -^ e(A - s\k4\) = 6(A - sKf(\A\ - 1)). (33) 



First we eliminate one variable K4 using the 6 function 5^(X4+Clearly for general values of | A| the 6 function doesn't scale 
q-i- K2- qi), then the boson-fermion interaction term can invariantly. However, when 
be written as 



52- 
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y, r d"K2 d°q^ d^qi 
^ J (2^(2^(2^ 



|A| = 1, 



(34) 



f{K2 + qi - qi, J'„)f{K2, (^{,^)b''(q3,aj'^,^)b{qi,(^^,^) 

■S f ^ I, /■ „ 0{A - \k4\) 



the 6 function is invariant since 0(A) = 6(A/s). For the cou- 
phng constants in condition |A| ^ 1 we follow R.Shankar's 
analysis with a soft cutoff 11221 : 



■9(A-m6iA-\k2\)6(A-m), 



where 



1^:4! = \K2 + qi 



■q3\ 



Kf 



(29) 



(30) 



e(A - M) * e 



-|A-,I/A 



(35) 



Then the rescaled 6 function in our boson-fermion interac- 
tion term becomes 



Functions 6(A - Iq^l), 6(A - \k2\) and 0(A - \qi\) transform 

back to their original forms in the same manner as the pure 0(A - sKf\\A\ - 1\) ~ e""^^"^'"^' = g-'VAllA|-i|g-(..-i)iVAl|Ahi|^ 
boson case. However, the function 0(A - 1^4 1) is quite dif- (36) 

ferent here since k4 is a function not just of ^2, ^3 and qi where A^a = Kf/A. Since A <k Kf, we have A^a » 1. 

We can see if |A| = 1, the soft cutoff transforms invariantly. 



otherwise, it doesn't matter since the couplings will be ex- 
ponentially suppressed in the limit A^a — > oo. Hence, after 
the scaling the boson-fermion interaction term can be writ- 
ten as 

* ^^' B' 2-i j {2n)D {2n)'^ {2ji)0 

' Hi fiT M-. ^ ^ -^ ^ -^ ^ ' 



for Q<\K\-Kf < Af/s^, 

1 Y^ r d'^K 



,/(/t,c.{).'<^-^""^> 



/'(4 + q\ - q\,Oj'l + cJl_ - <)/'(4, w;,{) 



h'\q'yJ:;)h'{q\,a>l)(d(K). 



Then we can identify 



Si 



s'-"82. 



that is, the scaling dimension of g2 is 



[^2] = 2 - D. 



(37) 



(38) 



(39) 



for Af/s^ < \K\ -Kf < A/. 

Then the partition function can be recast as 

X I D[0;,0>,i/r>,l/r>] 

.g-S>[0>,0>,iA>,iA>]-'5/[0i.0<.iA<.iA<.iA>.0>.iA>.iA>] 



(42) 



(43) 



We next construct an effective action by integration over the 
fast fields. To the one-loop order, one obtains 



e 



-S,ffli/t%,if'<,ilit,i//<] 



At tree level the scaling dimensions of coupling constants 
g\ and g2 are both 2-D. This agrees with the reference 11291 
for the pure boson interaction. Hence, in 2 dimension they 
are all marginal. 

3.2 One-loop analysis 

In order to carry out the first step of Wilsonian renormal- 
ization group transformation at one-loop level, we need to 
perform a functional integration over the high-momentum 
part in the action. For convenience we spht the fields into 
"slow modes" and "fast modes". 



— e 
■exp 
1 



-S<[(^<,0<,iA<,(i'<] 



- (5/[0* , <P<, ^<, ^<, 0* , 0>, iA>' 'A>])^ 
-(5/[0* , (^<, 0-<, iA<, 0>, 0>, 4f>, 'P>f')^ 



(44) 



and 



where 



(f,T) = 0<(f,T)-K^>(f,T) (40) 



i/r(f , t) = i/r< (f , t) + i//^ (x, t) (41) 



for < 1^1 < Ab/s, 

for Ab/s < \q\ < Ab, 



where (...) denotes the average over the fast fluctuations, 
we perform the integrals over the fast modes by evaluat- 
ing the appropriate Feynman diagrams contributing to the 
renormalization of the vertices of interest. The one-loop 
Feynman graphs contributing to the renormalization are shown 
in Fig.l. After the integration over the fast fields we per- 
form the scaling transformations q' — > sq, k' — > s^k, b'iq'-, <jj'„) - 
^(-D-2)/2^(^.^ w„) and f'(k\, u'„) -^ s-^f(ki, oj„), which bring 
the cutoff' A/ s back to A. To keep the action invariant under 
renormalization transformation one finds that the chemical 
potentials and the coupling constants scale according to the 
following relations up to one-loop order. 

Hb -^ s'^\^b-2gi ( -pr-^Nsieq - ^^b) 
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Figure 1: The Feynman graphs contributing to the renor- 
malization of (a) the boson chemical potential Hh, (b)the 
chemical potential of the low-energy modes of fermions 
6^f, (c) the boson-boson interaction, and (d) the boson- 
fermion interaction. Dashed lines denote the boson fields 
and solid lines denote the fermion fields. 
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(NB(eg-fih)+l)\, (48) 



where 



NB(£q - l^h) = 



1 



j/'(f,-/'/>) _ \ 



(49) 



and 



Npiypk - 6fif) 



1 



(50) 



are the Bose-Einstein and Fermi-Dirac distribution func- 
tions which result from the summation over the Matsubara 
frequencies wjj and (Xt„ and OP is the D-dimensional solid 
angle, setting s - e^ and Ai(0) = A/(0) = A, we obtain the 
Gell-Mann-Low equations: 
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(27t)D 
+Nf(vfA- diifU 
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26nf - g2 .- . r, NB(e\ - lJ.h), 
(Zn) 
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■[NB(eA - Mh) + U + 

2A 



4/3NB(eA - t^h) 

l+2NB(£A-^lh) 
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(27Tr 
■[NF(-VFA-6Mf)-l] 

H-A^F(Vf A - 5yU/)[A^f (Vf A - (JjU/) - 1] 
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(2 - D)g2 - 2gig2 
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■pNBieA - HbWBieA -jUi) + 1], 



d0 
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where eA - ■^. Eq.(55) ^ - -2y6 shows that for large ^ the 
temperature T{C) always flows to infinity for nonzero initial 
temperature. This means in the vicinity of the critical point 
the Bose distribution and Fermi distribution can be reduced 

as 
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(56) 



Nf(vf^ - Sfif) 
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/3{vfA - 6fif) + 2 2 
To absorb the factor e^^ in the Eq.(56) we redefine the scal- 
ing of the interaction coupling constants in Eq.(51)-(55) as 



gi(£) = e'-''-">'gi and g2i{) = e'''-">'g2. Then the Gell 
Mann-Low equations are approximated as: 
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We observe that the contributions of the fermion loops go 
to zero as ^ — > oo in above equations because of the factor 
e 2^^ and e '*'^. Hence, in the vicinity of the critical point we 
can ignore these contributions. If we redefine the chemical 
potentials and the coupling constants as 

Jib = t^bla, 
6fif = S/ib/a, 
8i =81/7, 
82 = g2/7. 



where a - — and y 



(2nfa^ 



equations can be further simplified as: 



(62) 
the the Gell-Mann-Low 
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Figure 2: Flow diagram of the running coupling constants 
g\ and g2 in 3 dimensional case. 
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{\i2-hr 



(65) 



(66) 



The first terms on the right-hand side of Eq.(63)- Eq.(66) 
are from the tree level scalings. Notice that the tree level 
scahngs of the coupling constants ^1 and ^2 go as 4 - D 
instead of 2-Z). This is because that near a classical critical 
point the quantum theory reduces to the classical theory. 
The same situation has been discussed by reference [28]. 

For instance, we consider the 3 dimensional case. The 
fixed points can be calculated as 



and 



Ol/6,%|i,|2) = (0,0,0,0), (67) 



(/ii,%|i,|2) = (^,0,^,0). (68) 

The first one is the trivial Gaussian fixed point and the sec- 
ond one is the Wilson-Fisher fixed point. Around the Wilson- 
Fisher fixed point, the running of the two coupling constants 
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are shown in the flow diagram Fig. 2. We can see that with in Eq.(2), we can combine the boson and fermion field as 



a small negative initial value the coupling constant g2 runs 
to negative infinity. This could indicate a boson-fermion 
pairing instability. 



4 Conclusion 



a doublet <1) 



-A 



which is called superfield. Then the 



action can be rewrite in terms of superfield as 

r rP t tf 

5[(D',(D]= I d'^x I drl^Hdr-- — V^-n)Cy 

J Jo ^ 2mh 



+ |(<1)1'(1))2 



(69) 



In fliis paper we investigated a weakly interacting boson- This action is invariant under supergroup U{l\l)f&\. We 

fermion mixture model by application of Wilson's renor- used renormalization group method to explore if there is a 

malization group analysis. This model includes one boson- supersymmetry fixed point where jji, = yU/ and gi = g2. The 

boson interaction coupling constant g , and one boson-fermion calculation of the Gell-Mann Low equations shows that our 

interaction coupling constant g2. At tree level RG analy- model doesn't exhibit such a fixed point. 

sis shows that the scaling dimensions of gi and g2 are both 

2 - D. That is, the two coupling constants are marginal 

in D = 2. Here one needs to notice that the derivation of 

the scaling dimension of g2 is under a condition of Eq.(34), 

without which we won't be able to compare the rescaled ^iong for useful discussions 

action with the original one in RG transformation. 

At one-loop level we derived the Gell-Mann-Low equa- 
tions and found that in these equations the contributions 
from the fermion loops went to zero exponentially as ^ ^ 
oo compared with the contributions of the boson loops. We 
simplify these Gell-Mann Low equations by ignoring the 
fermion loop contributions and solve for fixed points in 3 
dimensional case as an example. We found two fixed points. 
One is the trivial Gaussian fixed point and the other one is 
the Wilson-Fisher fixed point at which g2 vanishes. This 
implies that the boson-fermion interaction decouples at the 
critical temperature. We also drew the flow diagram of the 
coupling constants gi and g2 around the Wilson-Fisher fixed 
point. We observe that g2 goes to negative infinity with a 
small negative initial value. This can be a boson-fermion 
pairing instability. 

Supersymmetry is a symmetry that relates boson and 
fermion. It has been one of the most active research areas 
in the high energy phvsics OOl . Various researches were 
also conducted to find supersymmetry in condensed matter 
systems||7][8l. If we have jjh - fij - ^ and gi - g2 - g 
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